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Abstract 

We study the time evolution of the quantum fluctuations of the axion field 
for both the QCD axion as well as axions arising in the context of supergrav- 
ity and string theories. We explicitly keep track not only of the coherently 
oscillating zero momentum mode of the axion but also of the higher non-zero 
momentum modes using the full axion potential. The full axion potential 
makes possible two kinds of instabilities: spinodal instabilities and paramet- 
ric resonance instabilities. The presence of either of these instabilities can 
lead to a quasi-exponential increase in the occupation of non-zero momen- 
tum modes and the build-up of the quantum fluctuations of the axions. If 
either of these becomes a significant effect then axions would no longer be 
a suitable cold dark matter candidate. Our results confirm the conventional 
wisdom that these effects are not significant in the setting of an expanding 
FRW universe and hence axions are indeed cold dark matter candidates. 
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I. INTRODUCTION AND MOTIVATION 



Over the years there has been growing evidence that nonbaryonic, cold dark matter plays 
an important role in the structure and evolution of the universe ( see, e.g. Axions 
are among the most promising candidates for the non-baryonic cold dark matter of the 
universe The concept of the axion was originally introduced and developed to solve 

the strong CP problem of QCD in an appealing and phenomenologically acceptable way [|7| . 
In these models, axions are pseudo-Goldstone bosons of a U(l) symmetry. After the advent 
of supergravity and string theories, it became clear that particles with the properties of 
axions were in fact more generic, and that in these theories there could be additional axion 
fields which may play a significant role in the history of the universe (see e.g. ||). 

In this paper, we examine the time evolution of the quantum fluctuations of the axion 
field as it oscillates about the minimum of its potential. The standard picture is that the 
axion oscillates coherently in its potential; only the zero-momentum mode is important. 
The original papers presenting the axion as a dark matter candidate [|]-|6j considered simple 
estimates of instabilities in the axion field that could result in energy being pumped from 
the zero mode into higher-momentum modes. If such an effect were significant, then it is 
possible that the energy stored in the axion field would be largely converted to kinetic energy 
and subsequently redshifted away. According to the original estimates, these effects are not 
significant. 

However in recent years we have come to realize that there are two kinds of instabilities 
occurring in the time evolution of generic mode functions that have the potential of changing 
this situation. Thus, either spinodal instabilities [|9Hl2f or parametric resonance instabilities 
T3|-15|, if they last for a sufficiently long time, can lead to an explosive growth of quantum 



fluctuations through the exponential growth of non-zero momentum modes. Thus in the 
light of the recent understanding of the role of these instabilities in the growth of non-zero 
momentum modes, it is worthwhile to re-examine the role of the quantum fluctuations of 
the axion. 

In what follows we carefully and quantitatively study the time evolution of both the zero 
and non-zero momentum modes of the axion. We are thus able to ascertain the magnitude 
of the quantum fluctuations of the axion and compare it with the value of the coherently 
oscillating axion field. We do this both for the QCD axion as well as other axions that 
arise in the context of supergravity and string theories. We show that if the axions were 
born, lived and died in Minkowski space then there would in fact have been an explosive 
growth of quantum fluctuations resulting from the quasi-exponential growth of some non- 
zero momentum modes of the axion field. However, the energy density of axions is diluted by 
the expansion of the universe, which implies a decrease in the amplitude of oscillation of the 
coherently oscillating zero mode. Since it is this oscillation that drives the instabilities and 
the explosive exponential growth of the non-zero modes, it is clear that when the amplitude 
of the zero mode falls below some critical value, the instabilities will be shut off and there 
will be no further growth of the fluctuations. The issue thus becomes one of the initial 
amplitudes and timescales involved. 

In our analysis we will restrict our attention to the time evolution of the axion zero mode 
and axion fluctuations. We will neglect the couplings of the axions to other fields p!6|JT7[] . In 
principle, the coupling of the axions to say photons could result in a parametric resonance 
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production of photons. In fact, parametric resonance results from the transfer of energy 
from the zero mode of the axion in this case to modes of other fields. This can only occur if 
the axion zero mode oscillates with a sufficiently large amplitude for a sufficiently long time. 
The central issue here is whether in the cosmological context, the axion zero mode oscillates 
with a sufficiently large amplitude for a sufficiently long time to make parametric resonance 
possible. While in our quantitative analysis we will only keep track of the axion zero mode 
and axion fluctuations and neglect the couplings of axions to other fields we expect that 
if the axion zero mode is not able to excite the non-zero modes of the axion then it will 
also not be able to excite the modes of other fields. If in fact we did find that the non-zero 
momentum modes of axions were getting excited by parametric resonance then we would 
have to keep track of the modes of other fields. 

We now turn to a quantitative analysis of the problem to determine whether there is 
enough time to build up the fluctuations significantly. In the next section we will describe and 
layout the equations that determine the time evolutions of the axion field and its fluctuations. 

An important point is that we will keep the terms in the potential which are non-linear 
in the mean field <p/f a ] that is, we will keep the full expression cos (</>//„) instead of just 
the terms up to quadratic order in the mean field <p 2 . This is necessary since the initial 
value of the mean field can be of order one, and so neglecting the higher order terms in 
the potential is unwarranted. However, the main issue we wish to address in this paper is 
whether or not the fluctuations in the axion field ever become large due to the presence 
of instabilities in the evolution equations for the mode functions. For this purpose it is 
sufficient to start off with the initial fluctuations small, and then keep only the lowest order 
term in the fluctuations of the axion {rj 2 )/ f 2 , where r\ represents the fluctuations in the axion 
field about its mean value 0. If in following the time evolution we were to find that (rj 2 ) j f 2 
did become of order one, then we would be required to deal with the non-linear terms in the 
fluctuations in some way, such as the Hartree approximation ( see e.g. [0]). However, as we 
shall see, the fluctuations never become large in a cosmological setting, and so we will not be 
required to use the Hartree approximation. However, in order to arrive at this conclusion, 
it is nevertheless necessary to keep the non-linear terms in the mean field 4>/f a in the axion 
potential. Further, we also consider the explicit time dependence of the axion potential in 
a cosmological setting due to the temperature dependence of the instanton effects that give 
rise to the axion potential. These issues are discussed in greater detail in section 2. 

In section 3 we will study and analyze the solution to these equations for (i) the QCD 
axion and (ii) other axions arising from supergravity and string theories. In both cases we 
will compare the evolution in Minkowski space with the behavior in an expanding FRW 
universe to gain insight into the relative roles of instabilities and the expansion of the 
universe. Finally we will conclude by stating the implications of our results and place things 
in perspective. 



II. EVOLUTION EQUATIONS 



The derivation of the appropriate evolution equations [^,0] has been intensively studied 
during the past few years by a number of groups |T ^To]JTE| -pS[ . Here we will summarize the 
key formulae along the lines presented by Boyanovsky, de Vega and Holman |31 ] . 

In a spatially flat FRW cosmology, the metric is 
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ds 2 = dt 2 - a 2 (t)dx 2 



(2.1) 



The action and Lagrangian density for the axion field $ are given by 



S = J d 4 xC 
£ = a 3 {t) 

V($,T) = A 4 (T) 



!a2^ x l(V$(f,t)) 2 Tr/ _ , N 



1 — COS 



fa , 



(2.2) 
(2.3) 

(2.4) 



where f a is the axion decay constant and T is the temperature, and A(T) is related to / 
and the axion mass m a {T) via A 2 (T) = f a m a {T). The temperature dependence of m a and 
A is due to the temperature dependence of the instanton effects which give rise to the axion 
potential []3^ , |33| . 

For reasons explained in section 3, we will concentrate for the most part on the time and 
temperature dependence of the QCD axion mass, and in the case of axions arising in the 
context of supergravity and string theories, we will restrict our attention to the case of a 
time independent potential. 

For our purposes the temperature dependence of the QCD axion mass can be captured 



by the expression 32,|33|: 



m a (T) = 0.1 m, 



uO 



A, 



QCD 

T 



.3.7 



(2.5) 



where m a0 is the zero-temperature axion mass. It should of course be noted that this form 
of the above expression is only valid for T > Aqcd- 

We will be interested in two kinds of axions: the standard QCD axion for which Aqcd ~ 
200 MeV and f a ~ 10 12 GeV, and axions that can arise in the context of supergravity and 
string theories, for which we will take A ~ 10 16 GeV and f a ~ Mpi ~ 10 19 GeV (although the 
parameters are much less constrained in this case). Axions with these parameters have been 
considered earlier in the context of natural inflation ||34|| . We will now consider computing the 
magnitude of the fluctuations of the axion field compared to the amplitude of the coherently 
oscillating zero momentum mode of the axion for general values of A and f a . 

The canonical momentum conjugate to $ is 



U(x,t) = a 3 {t)$(x,t) , 



and the Hamiltonian becomes 



H(t) 



d x 



j n 2 

\2a 3 {t) 



+ 



a(t) 



(V$) 2 + a 3 (t)V($) 



(2.6) 



(2.7) 



In the Schrodinger representation (at an arbitrary fixed time t ), the canonical momentum 
is represented as 



II(£) = -% 



5 

8$(x) 



4 



Wave functionals obey the time dependent functional Schrodinger equation 



' dt 



HV[$,t] . 



(2.8) 



For the systems we'll be interested in it is convenient to work with a functional density 
matrix p with matrix elements in the Schrodinger representation p [<&(?), $(7); t}. Normalizing 
the density matrix such that Trp = 1, the "order parameter" is defined as 



i 



d 3 x(§(x, t)) 



n 



d 3 xTrp{t)${x) 



(2.9) 



where Q is the comoving volume, and the scale factors cancel between the numerator (in the 
integral) and the denominator. Note that we have used the fact that the field operator does 
not evolve in time in this picture. In this paper we will use the terms "order parameter" , 
"mean value of the field" and "zero momentum mode of the field" interchangeably to refer 
to quantity defined above. The evolution equations for the order parameter are 



d(f>{t) 

dt 
dir(t) 

dt 



a 3 (t)n 



d 3 x (U(x, t)) 



-— / d 3 xa 3 (t) 



a 3 (t)tt 
6V($) \ 
5<$>{x) I ' 



d 3 xTrp(t)U(x) 



At) 
a 3 (t) 



It is now convenient to write the field in the Schrodinger picture as 

$(£) = cf)(t)+r}(x,t) 
( V {x,t)) = 0. 



[2.10) 
(2.11) 



(2.12) 
(2.13) 



Expanding the right hand side of (|2.11|) in powers of rj(x, t) we find the effective equation of 
motion for the order parameter: 



v"W) ) r , u 

dt* 1 ~a(t) dt~^ V W»^—XT 



d 3 x(r) 2 (x,t)) + 0(ri 4 ) = . 



(2.14) 



where primes stand for derivatives with respect to 0. For our case, with V A ($) given by ( f2.4| ), 
we find 



d 2 (j)(t) a(t)d(f>(t) A 4 (T) . ( <j> 



dt 2 



+ 3 



a(t) dt 



+ 



f Sm \f 

J a \ J a 



1 



+ o( v 4 



o 



(2.15) 



It is legitimate to neglect the 0(r] 4 ) and higher terms in this equation as long as (rj 2n ) / f 2n 
remains small. We will assume that this is the case provided {t] 2 )/fa remains small. Of 
course, if we were to find (rf)/f 2 of order 1 at some time, then beyond that time it would 
not be legitimate to neglect the higher order terms, and we would need to keep track of 
them through some non-perturbative technique such as the Hartree approximation. The 
first issue to address is if and when (jf) / f 2 ever does become of order 1, and this we can 
do while neglecting the 0(^ 4 ) terms. To do this it is sufficient to consider only the lowest 
order terms in (f] 2 )/f 2 - 
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In order to follow the time evolution of the fluctuations (rj 2 ) it is convenient to introduce 
mode functions ipk{t) which obey a simple evolution equations: 



d*<Pk 
dt 2 



k 2 A 4 



(T) 



( 



fa 



) 



(v 2 )] 



(2.16) 



+ 3- 



a dt 



+ — + 

or 



f 2 

■J a 



COS 



1 - 



¥k = , 



where we have dropped the 0(f] A ) and higher terms. In terms of the functions <fk{t) the 
initial conditions are taken to be 



where Wk(t ) = [a~ 2 (t )k 2 + m 2 ^ 1 ^ 2 . This initial condition corresponds to taking an initial 
gaussian wave-packet for the field <3? with a width determined by the parameter mo which has 
the dimensions of mass. The smaller the parameter mo the more sharply peaked the initial 
distribution is around k — 0. Since we are interested in studying the process of building up 
the occupation of higher non-zero momentum modes starting off with essentially only the 
zero momentum mode we are interested in the case where mo is much less than m a . We 
have also varied m and checked that the rate of the build-up of fluctuations is insensitive 
to the choice of m®. This is actually so because when the non-zero momentum modes do 
grow they do so at an exponential rate and hence any differences in the initial occupation 
of modes quickly becomes insignificant. 

The equal time two-point function for the fluctuations which we have been denoting by 
(rj 2 ) can then be expressed as: 



At this point, it is worth discussing the momentum cutoff in the integral above. The 
potential for the axion is non-renormalizable and the theory we have described above should 
be considered only as a low energy effective field theory. Having said that we still have to 
address the issue of the cutoff in the momentum integral. First it should be noted that 
the instabilities in the mode functions can occur for momenta k less than or order of m a . 
Thus for instance for the QCD axion, if we choose the momentum cutoff k cu toff such that 
m a « k cuto ff << Aqcd then we will certainly capture any instabilities in the growth of the 
mode functions. In particular, we have taken the momentum cutoff ~ 10 2 m a and verified 
that our results are insensitive to the cutoff. Note that when the instabilities in the mode 
evolution equation actually kicks in (as happens for instance in Minkowski space) the mode 
functions which can undergo a quasi-exponential expansion are well within our momentum 
cutoff. 

Having described the formalism and arrived at the evolution equations for the axion 
fields of interest for us we are now in a position to study the solutions to these evolution 
equations. We turn to this in the following section. At this point we point out that there 
are some aspects of a similar problem specifically in the context of natural inflation that are 
currently being investigated [j35|l . 
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III. SOLUTIONS TO THE EVOLUTION EQUATIONS AND ANALYSIS. 



To analyze the solutions to the evolution equations we have written down in the previous 
section it is convenient to rescale variables into dimensionless ones in the following way. 

For a radiation dominated (RD) universe, the temperature dependence of the axion 
potential can be translated into an explicit time dependence by using |33] 



H{t) = 777 



1 

2t 
5 

3 9 



1/2 



CO 



rp2 

mpi 



(3.1) 
(3.2) 



where H(t) is the Hubble parameter as a function of time and g*(T) is the number of effective 
degrees of freedom at temperature T. The oscillations of the axion field in its potential 
begin approximately at a temperature T± such that 3H(Ti) = m a (Ti); for the QCD axion, 
T\ ~ 0.9 GeV. It turns out to be convenient to write out the evolution equations in terms 
of dimensionless quantities defined in units of m a (Ti) in the following way: 



£ = T' r = m a (T 1 )t, ( x 2 ) = ^jt, 7] k = fcV /m a (T!) , q= — . 

We further introduce the dimensionless coupling 

1 (m a {T x )\ 2 



(3.3) 



fa 



(3.4) 



In terms of these dimensionless variables, the rescaled equations are 



d 2 i{r) , o a(r)d£(r) 



dr 2 

72 



+ 3 



air 



d 



d" 

dr 2 a(r) dr 



dr 



+ 



+ (O.lK-'sinf 



(x 2 ) 



a 2 (r) 



(0.1)r^cos£ l-(x') 



r)k{r) = 



(X 2 ) = 9 / (fdq 



(3.5) 
(3.6) 
(3.7) 



We will study axions and their fluctuations in two different space-time settings. We 
do this in order to compare the behavior and gain insight into the physics of the build- 
up of fluctuations. Thus we will study axions in a radiation dominated (RD) expanding 
universe, and in a static Minkowski space-time. Both of these situations can be captured by 
parameterizing the time evolution of the scale factor as 



a r 



a r 



(3.8) 



We will take a = 1 and consider the two cases n = 1/2 and n = corresponding to an RD 
Universe and Minkowski space-time respectively. For each space-time setting we study two 
different kinds of axions: the QCD axion and axions that arise in the context of supergravity 
and string theories. These are parameterized by their different values of g. 
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Our goal here is to study the question of whether instabilities in the mode equations 
can lead to a significant build up of quantum fluctuations of axions. We are interested in 
two kinds of axions: the standard QCD axion and axions that can arise in the context of 
supergravity and string theories. We will concentrate most of our effort here in the study of 
QCD axions for 2 reasons: (i) the parameters for the QCD axion are more tightly constrained 
by phenomenological, astrophysical and cosmological considerations and (ii) the potential 
implications of the quantum fluctuations of QCD axions has a more pressing relevance for 
cosmological issues. Thus we will study the time evolution of the QCD axion taking into 
account both the expansion of the universe as well as the explicit time dependence of the 
axion potential in a cosmological context. We will see that the time dependence of the axion 
potential only helps in further suppressing the growth of axion fluctuations. Thus if the 
axion fluctuations are unable to grow in a time independent potential, they will not grow 
in a potential in which mass increases as a function of time. Hence for the axions arising 
in the context of supergravity and string theories we will only investigate the growth of 
axion fluctuations in a time independent potential. We will see that for these axions too the 
fluctuations in an expanding universe do not grow even in a time independent potential. 

These two different kinds of axions have very different values of A and f a which enter 
into the quantity g which arose in the way we rescaled variables into a dimensionless form. 
Thus for the QCD axion we have g = 3 x 10 -60 , and for string and supergravity axions 
we have g = 10~ 14 . In this latter case, as we have already explained, we will restrict our 
attention to the case of a time independent potential. The relevant equations are the same 
as those of eqs. (p.3 H3.7|) , but with m a (Ti) replaced by m a0 ; also, the factors of 0.1 r 3 7 are 
absent in eqs. (|3.5|-|3.7|). 

Having done all the above things, the solutions to the evolution equations can be dis- 
played in terms of plots of quantities as a function of time. We have organized these into 
four figures: Figure 1 has plots representative of the QCD axion with a time dependent 
mass in a RD universe (n = 1/2, g = 3 x 10~ 60 ); Figure 2 has plots representative of the 
QCD axion with a time independent mass in a RD universe (n = 1/2, g = 3 x 10 -60 ); Fig- 
ure 3 has plots representative of supergravity and string theory axions in a RD universe 
(n = 1/2, g = 10~ 14 ); Figure 4 has plots representative of supergravity and string theory 
axions in a Minkowski space-time (n — 0, g — 10~ 14 ). 

For each of these cases above we plot the time evolution of the axion field zero momentum 
mode (£); and the logarithm of the axion fluctuations (log[(x 2 )])- 

We wish to emphasize that the purpose of figures 2 and 4 is merely to provide a clearer 
picture of the different roles played by the time dependence of the axion mass and of a static 
Minkowski space versus an expanding FRW universe respectively. Thus figures 1 and 2 differ 
only in the fact that whereas figure 1 displays the evolution with the time dependent axion 
mass in the cosmological setting, figure 2 has the evolution with everything the same except 
that the axion mass is frozen at the value at which the axions would start oscillating. Thus, 
by comparing figures la and 2a it is clear that the effect of the increasing axion mass is to 
increase the frequency of oscillations and to cause a more rapid dying out of the amplitude 
of oscillations. Further, by comparing figures lb and 2b it is clear that the fluctuations 
die more rapidly in the presence of a time dependent axion mass. Finally, figures 3 and 4 
highlight the difference between the Minkowski space and the FRW universe behavior and 
importance of axion fluctuations. In particular, it shows that with a time independent axion 



S 



mass, axion fluctuations would quickly become significant in Minkowski space whereas they 
would not be significant in an expanding universe setting. It should also be pointed out 
that in fact the qualitative behavior for both the axions arising in the context of string 
and supergravity theories and the QCD axion are similar, however there are quantitative 
differences because the parameters are very different in the different cases. 

IV. CONCLUSIONS 

In this paper we examined the time evolution of the axion field and its quantum fluctu- 
ations. This was done by keeping track not only of the zero momentum mode of the axion 
field but also of the higher non-zero momentum modes of the axion. We studied different 
kinds of axions: the axion required to solve the strong CP problem in the context of QCD, 
as well as other axions that would arise in the context of supergravity and string theories. 
Further, we studied the dynamics in two different space-time settings: in Minkowski space 
as well as in an expanding radiation dominated FRW universe. 

This study was motivated by the recent progress in understanding the important role 
of phenomena such as spinodal instabilities and parametric resonance instabilities that can 
lead to a rapid build-up of fluctuations. Indeed, we found that in Minkowski space axion 
fluctuations rapidly build up in a quasi-exponential manner. In the expanding universe there 
are a few additional factors which we understand through our quantitative analysis. First, 
we note that because the axion is a pseudo-Goldstone boson of a compact U(l) symmetry, 
the potential for the axion is periodic and hence the displacement from the minimum is 
bounded. This limits the initial amplitude of coherent oscillations of the axion. Further, 
since the instabilities in the higher momentum modes are driven by the coherent oscillations 
of the zero momentum mode it is clear that if the amplitude of oscillations drop below some 
threshold the instabilities will not get a chance to really take off. The result of our analysis 
is that in fact in an expanding universe quantum fluctuations of the axion do not become 
significant. 

Indeed this result preserves a central piece of axion lore. If the higher non-zero mo- 
mentum modes had been significantly occupied and quantum fluctuations had not been 
negligible then it would not have been accurate to think of axions as oscillating coherently 
and axions would not have been a suitable cold dark matter candidate. 

ACKNOWLEDGMENTS 

We thank Rich Holman for helpful discussions. A.S. and M.S. were supported by the 
National Science Foundation grant PHY-91-16964. EWK was supported by the DOE and 
NASA under Grant NAG5-2788. 



9 



REFERENCES 



[1] Particle Physics and Cosmology: Dark Matter, M. Srednicki (ed.), Elselvier Science 
Pub., 1990. 

[2] The Early Universe, E.W. Kolb and M.S. Turner, Addison - Wesley, 1990. 
[3] G. Jungman, M. Kamionkowski and K. Griest, Phys.Rept.267: 195, 1996. 
[4] J. Preskill, M. Wise and F. Wilczek, Phys. Lett. B 120, 127 (1983). 
[5] L. Abbott and P. Sikivie, Phys. Lett. B 120, 133 (1983). 
[6] M. Dine and W. Fischler, Phys. Lett. B 120, 137 (1983). 

[7] R. Peccei and H. Quinn, Phys. Rev. Lett. 38, 1440 (1977) ; S. Weinberg, Phys. Rev. 
Lett. 40, 223 (1978) ; F. Wilczek, Phys. Rev. Lett. 40, 279 (1978) ; M. Dine, W. Fischler 
and M. Srednicki, Phys. Lett. B 104, 199 (1981). 
[8] T. Banks, M. Dine, "The cosmology of string theoretic axions", Nucl.Phys. B505 (1997) 
445-460; K. Choi, E. Chun and J. Kim, "Cosmological implications of radiatively gen- 
erated axion scale", Phys.Lett. B403 (1997) 209-217. 
[9] A. Guth and S.-Y. Pi, Phys. Rev. D 32, 1899 (1985). 

[10] E. Weinberg and A. Wu, Phys. Rev. D 36, 2474 (1987). 

[11] D. Boyanovsky and H. J. de Vega, Phys. Rev. D47, 2343 (1993). 

[12] D. Boyanovsky, D-S Lee and A. Singh, Phys. Rev. D 48, 800 (1993). 

[13] J. Traschen and R. Brandenberger, Phys Rev D42, 2491 (1990); 

Y. Shtanov, J. Traschen and R. Brandenberger, Phys. Rev. D51, 5438 (1995). 

[14] L. A. Kofman, A. Linde, and A. Starobinksy, Phys. Rev. Lett. 73 3195 (1994). 

[15] D. Boyanovsky, H. J. de Vega, R. Holman, D.-S. Lee and A. Singh, 
Phys. Rev. D51, 4419 (1995); 

D. Boyanovsky, M. D'Attanasio, H. J. de Vega, R. Holman and D. S. Lee, 

Phys. Rev. D52, 6805 (1995); 
[16] David B. Kaplan, Nucl. Phys. B 260, 215 (1985). 
[17] Mark Srednicki, Nucl. Phys. B 260, 689 (1985). 

[18] D. Boyanovsky, H.J. de Vega, R. Holman, J.F.J. Salgado, Phys. Rev. D54, 7570 (1996). 

[19] F. Cooper and E. Mottola, Mod. Phys. Lett. A2, 635 (1987); F. Cooper, S. Habib, Y. 
Kluger, E. Mottola, J. P. Paz, P. R. Anderson, Phys. Rev. D50, 2848 (1994); F. Cooper, 
S. Habib, Y. Kluger and E. Mottola, "Non-Equilibrium Dynamics of Symmetry Breaking 
in $ 4 Field Theory", Phys.Rev. D55 (1997) 6471-6503. 

[20] O. Eboli, R. Jackiw and S.-Y. Pi, Phys. Rev. D37, 3557 (1988); 
M.Samiullah, O. Eboli and S-Y. Pi, Phys. Rev. D44, 2335 (1991). 

[21] L. A. Kofman, A. D. Linde and A. A. Starobinsky, Phys. Rev. Lett. 76 1011 (1996). 

[22] E.W. Kolb, A.D. Linde and A. Riotto, Phys.Rev.Lett. 77 4290 (1996). 

[23] G. W. Anderson, A. D. Linde and A. Riotto, Phys. Rev. Lett. 77 3716 (1996). 

[24] E.W. Kolb and A. Riotto, Phys.Rev. D55 (1997) 3313-3317; A. Riotto and LI. Tkachev, 
Phys. Lett. B385, 57 (1996); S. Yu. Khlebnikov and LI. Tkachev, Phys. Rev. Lett. 77, 
219 (1996); S. Yu. Khlebnikov and LI. Tkachev, Phys.Lett. B390 (1997) 80-86. 

[25] D.T. Son, Phys. Rev. D54, 3745 (1996). 

[26] D.I. Kaiser, Phys. Rev D53, 1776 (1996). 

[27] H. Fujisaki, K. Kumekawa, M. Yamaguchi and M.Yoshimura, Phys. Rev. D53, 6805 
(1996); M. Yoshimura, Progr. Theor. Phys. 94, 873 (1995). 



10 



[28] D. Boyanovsky, D. Cormier, H. J. de Vega, R. Holman, A. Singh and M. Srednicki, 

Phys.Rev. D56 (1997) 1939-1957. 
[29] J. Schwinger, J. Math. Phys. 2, 407 (1961); P. M. Bakshi and K. T. Mahanthappa, J. 

Math. Phys. 4, 1 (1963); ibid, 12; L. V. Keldysh, Sov. Phys. JETP 20, 1018 (1965); A. 

Niemi and G. Semenoff, Ann. of Phys. (N.Y.) 152, 105 (1984); Nucl. Phys. B [FS10], 

181 (1984); E. Calzetta, Ann. of Phys. (N.Y.) 190, 32 (1989); R. D. Jordan, Phys. Rev. 

D33, 444 (1986); N. P. Landsman and C. G. van Weert, Phys. Rep. 145, 141 (1987); 

R. L. Kobes and K. L. Kowalski, Phys. Rev. D34, 513 (1986); R. L. Kobes, G. W. 

Semenoff and N. Weiss, Z. Phys. C29, 371 (1985). 
[30] For a thorough exposition of non-equilibrium methods in cosmology see, for example: 

E. Calzetta and B.-L. Hu, Phys. Rev. D35, 495 (1988); ibid D37, 2838 (1988); J. 

P. Paz, Phys. Rev. D41, 1054 (1990); ibid D42, 529 (1990); B.-L. Hu in Bannf/Cap 

Workshop on Thermal Field Theories: proceedings, edited by F.C. Khanna, R. Kobes, 

G. Kunstatter, H. Umezawa (World Scientific, Singapore, 1994), p. 309 and in the Pro- 
ceedings of the Second Paris Cosmology Colloquium, Observatoire de Paris, edited by 

H. J. de Vega and N. Sanchez (World Scientific, Singapore, 1995), p. Ill and references 
therein. 

[31] D. Boyanovsky, H. J. de Vega, and R. Holman, Phys. Rev. D49, 2769 (1994). 
[32] D. Gross, R. Pisarski and L.Yaffe, Rev. Mod. Phys. 53, 43 (1981). 
[33] For a discussion of this point and related cosmological issues please see section 10.3 of 
reference 0. 

[34] K. Freese, J. Frieman and A. Olinto, Phys. Rev. Lett. 65, 3233 (1990); F. Adams, J. R. 

Bond, K. Freese, J. Frieman and A. Olinto, Phys. Rev. D47, 426 (1993). 
[35] D. Cormier, R. Holman and S. P. Kumar, Work in progress: Private communication. 



11 



FIGURES 
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Fig. lb: QCD axion fluctuations with time dependent axion mass in a radiation-dominated universe. 
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Fig. 2a: QCD axion field with time independent axion mass in a radiation-dominated universe. 
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Fig. 2b: QCD axion fluctuations with time independent axion mass in a radiation-dominated universe. 




Fig. 3a. Time evolution of the supergravity and string theory axion zero momentum mode in a radiation- 
dominated universe. 
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Fig. 3b: Time evolution of the supergravity and string theory axion fluctuations in a radiation- 
dominated universe. 
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Fig. 4a: Time evolution of the supergravity and string theory axion zero momentum mode in Minkowski 
space. 
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Fig. 4b: Time evolution of the supergravity and string theory axion fluctuations in Minkowski space. 



